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MARTIAN TIDAL PRESSURE AND WIND FIELDS OBTAINED 
FROM THE MARINER 9 INFRARED SPECTROSCOPY EXPERIMENT 

by 

Joseph A. Pirraglia* 

Barney J. Conrath 

ABSTRACT 

Using temperature fields derived from the Mariner 9 infrared spectroscopy 
experiment, the Martian atmospheric tidal pressure and wind fields are calcu- 
lated. Temperature as a function of local time, latitude, and atmospheric pres- 
sure level is obtained by secular and longitudinal averaging of the data. The result- 
ing temperature field is approximated by a spherical harmonic expansion, retaining 
one symmetric and one asymmetric term each for wavenumber zero and wave- 
number one. Vertical averaging of the linearized momentum and continuity equations 
results in an inhomogeneous tidal equation for surface pressure fluctuations with 
the driving function related to the temperature field through the geopotential 
function and the hydrostatic equation. Solutions of the tidal equation show a di- 
urnal fractional pressure amplitude approximately equal to one half of the verti- 
cally averaged diurnal fractional temperature amplitude. These results indicate 

* NRC-NASA Resident Research Associate 

Preceding page blank 




that a diurnal pressure fluctuation of 6-7% existed during the planet-wide dust 
storm of 1971-72 as well as during its subsequent decay. The calculated tidal 
pressure fields, along with the temperature fields, yield tidal wind velocities of 
the order of 10 meters per second near the lower boundary, assumed to be fric- 
tion free in the model. Dynamic heating accounts for 3-4 °K of the diurnal 
amplitude of the vertically averaged temperature. 



MARTIAN TIDAL PRESSURE AND WIND FIELDS OBTAINED 
FROM THE MARINER 9 INFRARED SPECTROSCOPY EXPERIMENT 

INTRODUCTION 

Temperature profiles obtained from measurements in the 667 cm' 1 C0 2 
band by the infr ared interferometer spectrometer experiment on Mariner 9 have 
been used to construct semi-global temperature fields for the Martian atmosphere 
in the form of longitudinal averaged temperatures at various pressure levels 
given as functions of latitude and local time. During the great dust storm of 
1971, a strong local time dependence of the temperature throughout the full 
depth of the lower atmosphere was observed. Using a dynamical model which 
assumed a constant surface pressure, the temperatures fields have been em- 
ployed to estimate the Martian winds which existed during the dust storm (Hanel, 
et al, 1972). However, the temperature data suggest an appreciable tidal motion 
forced by diurnal heating. The object of the present paper is to investigate the 
tidal surface pressure variations and associated wind fields again using the ob- 
served temperatures as input data. In addition to the purely tidal motion, the 
investigation includes local time averaged fields, i.e., those fields which are 
symmetric with respect to the polar axis. 

SURFACE PRESSURE TIDAL EQUATION 

The set of equations which describes the motion of the atmosphere is sim- 
plified when the temperature, rather than the heat input, is used as the known 
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driving mechanism. The heating required to be consistent with the observed 
temperatures may then be calculated. Using the perfect gas law and the hydro- 
static approximation, the equations of motion in a- coordinates (Phillips, 1957) 
can be written 


3 V 1 I • 3a RT r~7 V 7 >k C I 

— + v • Vv + a — + — Vp +V$+fkx 

3t 3a p s 


v = F 


P ^=- Tv. 


3p 

(p s v)da- (a- 1) — 


$ - - R 



T 

— da' + 3> 
a' 


The following notation is employed: 

0, colatitude; 

0, longitude measured from west to east; 
fi, rotational angular velocity; 
t , time; 

R, C0 2 gas constant; 
p s , surface pressure; 
p , pressure; 
a, pressure ratio, p/p s ; 
k , unit vector in radial direction; 

V, horizontal gradient operator on constant a surface; 
v, horizontal velocity on a -surface; 


( 1 ) 

(2) 

(3) 
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T, temperature; 

a, vertical velocity in a- coordinates; 

F , frictional forces; 
f , coriolis parameter, 2 Q sin d ; 

$ s , the geopotential on surfaces = 1. 

Assume that the temperature may be expressed in terms of T q , a tempera- 
ture field associated with a static atmosphere, and T 1 , a temperature associated 
with atmospheric motion, in the form 


where e is a small quantity. Expanding the variables v,.<j, p s and <I> in terms 
of e about the static case, designating the static terms by subscript 0 and first 
order terms by subscript 1, and using the boundary condition cr = o at a = 0 
and cr = 1, the first order equations are 


Bv 


RT„ 


— 1 + fkx Vj +— °Vp sl + 1^ - 


T oPs\ V P S 0 


+ VOj = P’ 


sO 


sO / SO 


3 Psi 


Bt ' " s0 J 1 
o 


J' 


+ P.n V-v.da + Vp, 


sO 




dcr = 0 


(4) 

(5) 



( 6 ) 


Neglecting orographic effects, p^ is constant on the surface a = 1 and the 
fourth term of (4) and the third term of (5) can be dropped. For a static case 
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in which p sQ is not constant, the terms could be dropped and later used to assess 
the validity of neglecting orographic effects. Only the first order equations will 
be considered. 

Averaging (4) and (6) over cr from 0 to 1 and denoting a -averaged quantities 
by an over-bar, the equations of motion for the vertically averaged velocity and 
surface pressure are 


Bv, RT n 


+ fk x v + Vp + VO _ F 

Bt 1 p sl 1 

*s0 

(?) 

Bp 

it + p -» v ' *. • a 

(8) 

Oj = RT r 

(9) 


The friction term F depends upon the internal velocity shearing stress and 
the shearing stress at the boundary. The concern here is with the wind regime 
outside the surface boundary layer and with the surface pressure variations. 
Since the motion is forced, the bulk of the flow and the resulting surface pres- 
sure should not be greatly affected by the surface drag except possibly for the 
axially symmetric flow in which friction causes a small meridional component 
that allows a latitudinal redistribution of the atmosphere. With these observa- 
tions in mind and with the intent of simplifying the equations, the friction term 
is replaced by a linear damping term proportional to the velocity, F = Sv , where 
S is a damping coefficient to be specified. With this approximation the set of 
equations 7, 8 and 9 are complete and formally may be solved for 
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the velocity Vj and the surface pressure p sl , assuming that the static pressure 
and temperature, p sQ and T 0 , and the motion inducing temperature field are 
known. From this point on the subscript 1 denoting the first order terms will 
be dropped. 

Letting the first order temperature field be represented by a superposition 
of fields of the form 

T = T wk (0, CT ) ei( “ l + k4>) (10) 


solutions of equations 7 and 8 for v and p s of the same form will be sought. 
With At = cos 6 and the definitions, 


cok 


p T 

co k OM 

Pso T q 


( 11 ) 


4> 


co k 


co k y 

1 n 


(12) 


r - i r = —_i— 
^ - W + W ~ 2 n ~ 1 2fi 


(13) 


CO 


- co - i 8 


(14) 


equations (7) and (8) can be reduced to 


3 

"1 - At 2 

i 

k 2 k (£ 2 + At 2 ) 

3 At 

•Q - At 2 _ 

£ 2 -m 2 

_1 - At 2 ^ (£ 2 - M 2 )_ 


where a is the radius of the planet. 


4a 2 fi 2 

+ — = — 
RT n 


— ; G 


co k 


4a 2 fl 2 qj — » 

RT a? " k 

(15) 
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Equation 15 is an inhomogeneous Laplace Tidal equation with a linear damp- 
ing term. The homogeneous equation, without damping, has for its solutions the 
Hough functions which are the ^eigenfunctions associated with eigenvalues gen- 
erally expressed in terms of "equivalent depths" (Chapman and Lindzen, 1970). 
This problem has been discussed extensively in the literature (see for example 
Longuet- Higgins, 1967). The temperature field T could be expanded in Hough 
function and the inhomogeneous equation solved in terms of the eigenfunctions. 
However, in this case the eigenvalues and eigenfunctions would be complex due 
to the damping term. Rather than expanding in Hough functions which are infinite 
series of Legendre functions, the temperature field and resultant pressure and 
wind fields will be expanded in terms of Legendre functions. For the problem 
under consideration, the Legendre function solution converges reasonably fast. 
The only difficulty is that in solving the momentum equation for the wind veloci- 
ties there is an apparent singularity when damping is not included. A zero in 
the Hough function formulation cancels the pole, while in the Legendre function 
formulation each Legendre mode does not exactly cancel the pole but instead the 
sum of the modes is required. Since the solution in terms of Legendre functions 
converges rapidly and because of the damping term, the difficulty is a minor one. 

For the wave number k the temperature component T is expanded in a 
series of Legendre functions (m), 

m 

co 

t -k =2Z /3kmCc7)P '" (M) (16) 

m- k 
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where the vertical structure is contained in the coefficients /3, (o-). Using 

km 

equation (16) the geopotential function becomes 


m = k 


( 17 ) 


where /3* m = /3 km /T If the surface pressure is assumed to be of the form 



(18) 


then from (11) and (12) 



m= k m= k 


Now substituting (17) and (19) into the tidal equation (15) the recursion relation 
for the complex amplitudes is obtained; 


m 

i - k 

m-k-l_* i (m-l) 2 (m + k) 

(2m - 1) 

" , . k" 

m (m 

r\ 'X km - 2 ^ n 

^ m 2 (2m - 1) 


“kn 


m(m + 1) - — tv? 
r2 £_ 0 (£_ 

m 2 (m + l) 2 4a 2 n 2 w 

m + k + 1 


km 


( 20 ) 


(2m + 3) 


(m+ l)(m+ 2) -y-j 


(m + 2) 2 (m- k + 1) * , m + k + 2 


(m+ l) 2 (2m- 1) 


km 


2m + 5 C ' km + 2 


RT, 


0 w £ 


4a 2 fi 


2 n 2 co 


km * 


Equation (20) represents an infinite set of equations for the coefficients 
If the solution as represented by (18) converges then the set of equations can be 


a, • 

km 
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truncated to give a solution to the desired accuracy. The set of equations can be 


simplified somewhat by the introduction of auxiliary functions , 


(m + l) 2 (m - k) 
m 2 (2m - 1) 


a 


km — 1 


m + k + 1 
2m + 3 ° km+1 


m + 1 
m 



(21) 


Defining the following expressions; 




_ (m + l) 2 (m - k) 
(2m - 1) m 2 

(22) 

. _ m + k + 1 

Lkm ~ 2m + 3 

(23) 

, m + 1 

U m „ 

1 k 

n 2 £ 

RT 

(24) 

m + 1 lk 
m m 2 

„ o 

- (m + l) 2 

4a 2 f) 2 

(25) 


and substituting into (20) and (21) the equations for the coefficients are, 


^km'Atm- 1 


- o + 1 > 2 

a km- 1 “ ^km^km + ^km “km + 1 


RT 0 oj ' 75* 
4a 2 f? “ km 

(26) 

= o. 

(27) 


Equations (26) and (27) form two sets of equations, one set for odd and one 
set for even values of k-m. The even k-m modes are symmetric about the equa- 
tor and odd k-m modes are asymmetric. In practice, (26) and (27) can be trun- 
cated and expressed in matrix form to solve for the finite set of a km which 
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adequately represent the solution to the tidal equation. The approach used here 
closely parallels that used in the eigenfunction solution of the homogeneous tidal 
equation (Chapman and Lindzen, 1970). 

Having solved for the a* m the pressure fluctuations are obtained from the 
assumed form in (10) and from (18). 


00 


= Po e 


i(<i)t+k0) 


L 


km 


p>) 


(28) 


m = k 

with summation over the finite set of a* P k in practice. For multiple values of 

km m 

frequency co and wave number k the solution is represented by a linear super- 
position of terms of the form in (28). 


POLAR SYMMETRIC FIELDS 

The polar symmetric temperature field and the resultant pressure and wind 
fields can be regarded as zero wave number tidal theory. Because the time 
scale of any changes of these fields is much larger than the time scale of one 
day, the zero frequency limit of tidal theory can be used. 

If k is made equal to zero in (20) where damping is implicit and the zero 
frequency limit is taken, the result is 


* — ft* 

^Om Pom 


(29) 


independent of the magnitude of the damping factor. On the other hand, taking 
the limit with S equal to zero (no damping) the result is, 
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(2m- l)(2m - 3) 


a 0m- 2 + 


RT„ 


4s 2q 2 (2m+3)(2m-l) m(m + 1) 


Om 


l 


RX 


O o* 


(2m + 3) (2m + 5) 


Om + 2 


4a 2 fi 2 


yS* . 


(30) 


Thus, the limit of S and approaching zero is not uniform. In the harmonic case 
with finite, the solutions are continuous functions of the damping parameter. 

Considering the damped case as an initial value problem and applying the 
final value theorem, the result is identical to (29). The final value theorem does 
not apply to the undamped case. Assume that the initial velocity and pressure 
distributions are v(0) and p (0) in (7) and (8). With the temperature applied as 
the step function (T t = 0, t < 0; T = T (9 , o), t > 0} and taking the Laplace 
transform of (7) and (8) using (9) and the linear damping term, the following 
set of equations is obtained, 


sVj - Vj(0) + fkx Yj 


0 o R VTV 

+ Vp + S V = - _VT 

p „ 1 1 s 

^sO 


(31) 


^sl -Psl(°) 


p v- 

^s0 


Vj = 0, 


(32) 


where s is the complex frequency and the tilde over a function represents its 
Laplace transform. 

Equations (31) and (32) can be reduced to an equation analogous to the tidal 
equation (15); 
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_3_ / 1 - M 2 ^P* \ + s 4a 2 n 2 ^ _ JL_ 3 / 1 - /J 3$*\ 

^ U 2 - M 2 ^ / s + 8 RT 0 P s ^2 _ fj2 ^ J 


s 4a 2 O 2 1 /m 1 s + S 

— p (0) + 

s + § RT q L s s (s + S) 2 + f 2 


( v . <0) -rb k><v ‘ (0) )_ 


(33) 


where p* is the Laplace transform of p sl /p sQ , 4>* = T/T q , and i 2 - - (s + S ) 2 /4f) 2 » 
As was done previously 4** and fs” * are expanded in Legendre functions, i.e., 


=]2 ^ P » 

m= 1 



(34) 

(35) 


Substituting (34) and (35) into (33), the recursion relationship for the coefficients, 
a* , of the expansion of the Laplace transform of the pressure ratio is obtained; 


1 . _ s + S RT o 1 _ £ 2 ~ 1 a* 

(2m - l)(2m - 3) m ” 2 _ s 4a 2 Q 2 (2m+3)(2m-l) m(m + 1)J 

, 1 a* = - S + 8 ~j3* + (initial value terms) 

(2m+3)(2m+s) m+2 s 2 4a 2 n 2 


(36) 


Assume for the moment that the initial values are zero. Let the coefficients 


a* and B* be written as column matrices of vectors in the form 

m m 



ll 



Then (36) can be written in matrix form as 


A a = - 


s + 8 


RT„ 


s 2 4 a 2 Q 2 


P. 


(39) 


where A represents the matrix with diagonal and off-diagonal terms defined by 


W, 


s + 8 ^0 


£ 2 -l 


s 4 a 2 q 2 m(m + 1) (2m + 3) (2m - 1) 


(40) 


[A] 


1 


m, m+ 1 


(2m + 3) (2m + 5) 


(41) 


[A] '= 1 ( 42 > 

m ' m_1 (2m - l)(2m - 3) 

All other elements of A are zero. Using the final value theorem, (Gardner and 
Barnes, 1942) the time asymptotic solution for the inverse transform of the 
coefficients a* is. 


a= lim sa= lim 
s -0 s —0 



( 43 ) 


It is assumed that the infinite matrix A can be truncated at a point which will 
adequately describe the transforma and the inverse transform a . 



Taking the limits s - 0, the diagonal terms of A behave as (8 /s) RT 0 / 4a 2 Q 2 . 
The diagonal terms of A -1 behave then as (S/s) -1 4a 2 fi 2 /RT 0 and the off-diag- 
onal terms as (8 /sp (4a 2 R 2 /RT 0 / with j > 1. Using these limits 

= (44) 

s-0 s 4a 2 O 2 

where I is the identity matrix. Therefore (43) reduces to the same results as 
expressed in (29). 

If the initial values of p sl andp are not zero they will contribute terms to 
a of the form, 


lim A 1 
s-0 


a P 


s + S 


a 2 + 


(s + S) 2 + f 2 (s + S) 2 + f 2 


(45) 


where a lt a 2 and a 3 depend on the initial values but are independent of s. In 
the limit these terms go to zero, and the final state is independent of the 
initial state. 

To be physically realistic some dissipation should be present, and with 
the realization that the method by which the dissipation was introduced here may 
be oversimplified, the solutions to the polar symmetric equation will be assumed 
to be expressed by (29). Using this result along with the solutions to (26) and 
(27), the surface pressure can be obtained from the zero order quantities and the 
imposed first order temperature field. 
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VELOCITY FIELD 


The calculated surface pressure and the imposed temperature field provide 
sufficient information to calculate the horizontal wind field as a function of lati- 
tude, local time and altitude. With Vp sQ equal to zero and the friction term F 
replaced by the linear damping term-Sv 1 , the momentum equation (4) becomes 


— - + f k x v, + 8v = -Vp + V<j> (46) 

Bt 1 P. sl 1 

For harmonic fields of the form e V" t + k ^), using (10) and (16) and the definition 
/3* m = /3 km /RT 0 , (46) written in the form of the harmonic components \ ak of 
Vj is 



0 and <D are the unit vectors in the directions of increasing colatitude and east 
longitude. Taking v e wk as the 0 directed wind component, v^ wk as the 4> 
directed component and using the definition (13), the expressions for the wind 
field from (47) are 
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1 


RT, 


V, ■'■’k 


2Qa [ (fj 2 - + i 2 { ) 2 + 44 r 2 cf] /rr 




O 2 - C 2 + C 2 ) 


£ r kP k -^(1 -AX 2 ) 


£ (<- * /O x 

m=k 


■ - 2 £ r C 2 


dP 1 


d/j. 


kP v 


+ 2 H r £i 


C r kP^ (1 -M 2 ) 


dP k 


d/j. 


+ iC^ 2 -t? + if)kC i P“k<48) 


ET, 


^0 ,<wk 


2afl[( A i 2 -^ + £ l ?) 2 +4££ 2 i ] 


i= 2Z <«*„ + >< 


f dp k 


dP k 

(-&.»-{?♦£?) < 1 -^) 1 , d ; 

- 2 ^ 



dP k 

■2i^ 2 (W 2 )-^+i(M 2 -C? + tf) 


- M kP k + ^(1-M 2 )- 


dP 1 


dfj. 


(49) 


Expressions (48) and (49) are superpositions of products of ^-dependent functions 
and fj. (or latitudinal)-dependent functions. The ^-dependent functions are de- 
fined as being zonal or meridional modes. The zonal modes are, using a prime 
to indicate the derivitive of the Legendre function 


' r , kn 


/U? [(ax 2 -C 2 + C 2 ) 2 + 4 ^C 2 ] 


0* 2 -£ r 2 + £?) [^ r kP^-M(l-M 2 )P^'] 


-2^ r C 2 kP k| 


(50) 
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O 2 - ^ + ^ )k ^P 


* i , km 


A-fj . 2 [(jj? - i) + l 2 ) 2 + 4 ^ ^ 2 ] 

+ 2C r EiU r kP; -M(l -M 2 )P^']}- (51) 


The meridional modes are, 


M 


-(/x 2 -^ r 2 + ^ 2 xi-M 2 ) £ a p; 


r , km 


/l - M 2 Um 2 - £ 2 + £ 2 ) 2 + 4 £ 2 £ 2 ] 


-2C r Ci [-MkP^ + ^(lV)P^']| < 52 ) 


M. 


j- 2£ r £ 2 (1 ~^ 2 )Pm + Cm 2 - ^ 2 + ^ 2 ) X 


i , km 


/l-/x 2 [O 2 - £ 2 + £ 2 ) 2 + 4 £ 2 £ 2 




(53) 


The i subscripted zonal modes and the r subscripted meridional modes are due 
to the damping term and would be identically zero if there were no damping; the 
remaining modes would be simplified accordingly. 

The cr-dependent functions or mode amplitudes are defined as 


r , km 


RT 0 

2 aO 


[R e ( a * m ) + /3J 


(54) 


A. 


RT 


i - km 2 an 


° l«n( 0 - 


(55) 


16 



The real and imaginary parts of a* m are due to the complex coefficients in (26) 
and (27); the complex coefficients are again due to the damping and without 
damping A; kra is identically zero. Using the modes and mode amplitudes, (48) 
and (49) can be rewritten in condensed notation 


cf> , CO k 


-/ [A . Z . - A. . Z. . + i(A . Z. . + A. . Z . )] • 

/ r , km r , km i , km i t km v r , km i , km i , km r , km ' 

m= k 


(56) 


6, cok 


' y ’ "A r km M rikm A. (km M. jkm + i(A r kn) M i km + \ , km M r , km )] 


(57) 


m = k 


In the polar symmetric case where k = 0 and £ r = 0, the modes are: 


7 — — ^ ^ ^ pk • 

urn ~ 




**» = - 


/l - u 2 


r ok » 


^ 2 + i* 


£.P k/ . 

^ l m 


(58) 


(59) 


If there is no damping £. = 0, and the meridional mode is identically zero. The 
existence of the meridional polar symmetric wind depends upon the existence of 
damping, with the magnitude of the velocity depending on the damping parameter. 
The mode amplitude is from (29) 


RT 0 


( 60 ) 


and the polar symmetric winds can be written 
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V e , Om A)m ^Om 


( 61 ) 


6, Om 


V z ' 


im Om 


( 62 ) 


The total wind field can be expressed using the expressions defined above, 


F. (Ar .km Z r,km “ A i,k m Z i , km> COS C" + k <^ + 

co k m = k 

~j2 (Ar ’ kmZi 'k™ + A i.km Z r, km ) sin(cot + k 0 + 0 o>a)k ) +J2 A)nA 

m 25 k m = 1 


( 63 ) 


= FF. F, (k r.laA.km “ A i,km M i,km) + ^ + ^0,^ 

a) k m= k 

- J", (A r,km M i.km + A i,km M r,km) S in(crf + k<£ + „ k ) + A) .A 

m= k m= 1 


( 64 ) 


The angle 4 > is determined from local time of the temperature maxi- 
mum for a particular set of values of co and k. The co and k summations are over 
the combinations of frequency and wavenumber used to describe the temperature 
field. 


HEATING AND VERTICAL VELOCITIES 

The atmospheric heating consistent with the imposed temperature field can 
be calculated from the tidal winds and pressures obtained in the previous sec- 
tions. The heat equation in cr- coordinates can be written 
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( 65 ) 


3T • 3T RT / • Bp s „ 

Q = c — + c V • VT + C a — P.cr + a r—~ + cry • vp 

. p B t p p Bct o-p l s B t E 


and the vertical velocity a is determined from 


.-i-f 

Pc J, 


V • (p v) da - 


(a - 1) 3 P S 


Bt 


(66) 


Expanding in terms of the parameter e, the first order equations from (65) and 
(66) are 


Qi = 




^0 ^Psl 

P S 0 Bt 


(67) 


r<r 


a i=~ 


v 


Vj da - 


(a - 1) Bp s i 

P s o 


( 68 ) 


From this point on the subscript 1 will be dropped, and the first order quantities 
will have no subscript. 

In terms of the mode and amplitude functions defined in the previous sec- 
tion, the vertical velocity due to the polar symmetric fields are, from (62) and (68) 






(69) 


From (18), (56), (57), and (68) the vertical velocity associated with the diurnal 
motion is, to first order 
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Using the vertical velocity calculated from (69) and (70) the vertically averaged 
heating is calculated using (67). For the polar symmetric components the time 
derivitives are zero and the heating in terms of the mode functions is 



The heating associated with a diurnal wave of wave number k and frequency « is 



where p k is obtained from (18). 

In (72) both a wk and p s wk are complex. Q„ k is then the sum of real, 

Q r wk » and imaginary parts, Q. ak , which when combined with the exponential 
frequency and wavenumber term gives the expression for the total heating. 
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Q [Qr i<tfk cos (cut + kcp + <P 0iQ)k ) - Q i>fflk sin (wt + k 0 

oj k 


+ ^O.oil,)] + Qqi + Q 


02 


( 73 ) 


The above expression is the heating in ergs gm _1 sec -1 averaged over a ; 
multiplying this expression by the quotient of the surface pressure and the 
gravitational acceleration gives the heating in terms of ergs/cm 2 sec. 



(74) 


Using the hydrostatic equation in a -coordinates, B z/oo = - RT/gcr, the vertical 
velocity in cm/sec is obtained as a function of cr and the rate of change of 
temperature; 


dz 

dt 


To first order in e , 


dz 

dt 


RT R r A dT 

gcr g \cr dt 



cr 


d O'. 


RTo.J T/l 3T 

gO- g Jj \cr Bt 



The required value of 6- can be calculated from (69) and (70) 


(75) 


(76) 


= °01 +CT 02 + 



^r,o>k COS ( Wt + k( ^ + ^O.eok) 

-o- if « k sin(o>t +k0 + 0 o>(ak )] 


OJ k 


(77) 
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and 3T/3t from (10) can be written 


Sin + k< ^ + ^O.a,k)- ( 78 ) 

k 

APPLICATIONS TO MEASURED TEMPERATURE FIELDS 

The derivations of the previous sections will now be applied to the tempera- 
ture measurements from the Mariner 9 infrared spectroscopy experiment. In 
the results to be presented here, the mean surface pressure was taken to be 5 mb, 
and the temperature distribution was assumed to be of the form 

T (o* , fj., (p, t) = T 0 (cr) + 0 Ql (cr) P jOi) + /3 02 (cr) P 2 (fi) 

+ [yS u (o-)P} 00 +/? 12 (a)p2 (/U)] cos (fit + 0 + 4 > 0 ) (79) 

The height dependent coefficients were chosen by fitting the model to the tem- 
perature field measured during dust storm conditions. The resulting temperature 
contours at approximately one scale height are shown in Figure 1 and can be 
compared with the measured temperature field given by Hanel, et al. (1972). 

Since the data cover only the local time interval between approximately 5:00 and 
19:00 Martian local time, the model provides an interpolation for times outside 
this interval as well as for other data void regions in the northern hemisphere. 
Values of /3* m and /3* m (cr) are given in Table 1; these parameters are related 
to the coefficients in (79) by /3* m = /3 km /RT 0 , and /3* m is the aver aged value 

of PL (-). 
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The values of /3* m are used in (26) and (27) to solve for the coefficients a* m 
of the Legendre function expansion of the pressure ratio defined by (18). The 
values of a through a 16 , with and without damping, are listed in Table 2; only 
the first six modes are used since they converge rapidly. Using these a* m > the 
5 mb average surface pressure and the harmonic form suggested by (10), the 
diurnal component of the pressure field is obtained (Figure 2). With damping 
the lowest pressure occurs between the equator and the latitude of the diurnal 
temperature maximum and lags the temperature maximum in local time. 

Without damping the results would differ somewhat. The imaginary parts, 
which cause a phase shift in local time, and the real parts with even m, which 
are the equatorially asymmetric modes, are all identically zero. Therefore, 
the pressure low and vertically averaged temperature maximum would occur 
at the same local time and the pressure minimum would be approximately on the 
equator. 

The polar symmetric pressure component is determined by /3* m from (29). 
Adding the polar symmetric and the diurnal components, the surface pressure 
shown in Figure 4 is obtained. 

The modes of the horizontal wind field are calculated using (50), (51), (52), 
(53), (58), and (59). Examples of low order axially symmetric modes are shown 
in Figure 4 and higher order diurnal modes in Figures 5 and 6. In all cases the 
dotted curve would be zero if there were no damping and there would be singu- 
larities in the solid curves at ±30° latitude in Figures 5 and 6. 
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The amplitude functions defined in (54), (55), and (60) are shown in Figures 
7 and 8. From (54), (60) and table 2 the amplitude of modes 01, 02 and the real 
parts of modes 11 and 12 are linear functions of the difference in and all, 
or a large fraction of, its vertical average This situation occurs because 

these modes are driven by the temperature field directly, manifested in the geo- 
potential term, and indirectly through the surface pressure. All the other ve- 
locity modes are driven only by the surface pressure modes. Due to the behavior 
of these amplitude functions the components of the wind field associated with 
them change sign in the middle range of a- values. Since these are the dominant 
modes, the horizontal wind pattern reverses direction in going from the surface 
to high altitudes. 

From (63) and (64) the wind fields are calculated using the amplitude and 
mode functions. The near surface winds are shown in Figure 9, the one scale 
height winds in Figure 10 and the high altitude winds in Figure 11. In comparing 
Figures 9 and 11 the wind reversal is obvious. Near the surface the average 
zonal flow is to the west in the northern winter hemisphere and to the east, but 
diminished, in the southern summer hemisphere. The average surface meri- 
dional flow is from north to south in both hemispheres but stronger in the 
northern hemisphere. At higher altitudes these averaged flows are reversed. 

Vertical velocities are calculated using (76) with the results shown in 
Figure 12. The vertical velocity contours are similar but shifted 90° from the 
surface pressure contours. The dependence on the surface pressure is apparent 
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from the last term of equation (66). The dependence on the divergence of the 
velocity is more complicated and somewhat weaker . The effect of the polar 
symmetric fields is smaller than the diurnal components and is not readily ap- 
parent except near the poles. The polar symmetric vertical flow is downward 
from the north pole to within 10° north of the equator then upward to approxi- 
mately 75° south and then downward again to the south pole. This behavior is 
due to a combination of two symmetric pole to equator meridional cells and a 
pole to pole cell. 

The total atmospheric heating per vertically averaged column as calculated 
from (73) and (74) is shown in Fig. 13. This represents the total heating neces- 
sary to maintain the observed temperature field and calculated velocity field. It 
is of interest to examine separately, the dynamic contributions to the heating. 

The contribution due to vertical motion which is the second term in the right 
hand side of (67) is given in Figure 14, while Figure 15 shows the contribution 
due to the time rate of change of pressure or third term on the right hand side 
of (67). The dynamic contributions are both small compared to the total heating. 

In order to account for the required atmospheric heating at sub-solar latitudes 
(20°S) and local noon, the energy deposited in the atmosphere would amount to 
approximately 20% of the solar flux incident normally on the top of an atmospheric 
column. In order to achieve adequate night time cooling, an effective atmospheric 
emissivity of 0.53 would be required, assuming radiation to space (with an ef- 
fective temperature of 200K) is the principle energy loss mechanism. 
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SUMMARY 


The major dynamic effect produced by the observed diurnal temperature 
variation is the generation of diurnal pressure fluctuations (12% peak-to-peak) 
which are large by terrestrial standards. The resulting diurnal winds near the 
surface have amplitudes no greater than about 20 m/sec. Therefore, these 
winds alone are not capable of raising dust (Hess, 1972; Sagan and Pollack, 1967) 
and could not sustain the dust storm alone. However, by augmenting the polar 
symmetric wind fields or orographic wind fields, the diurnal winds could con- 
tribute to the raising of appreciable amounts of dust. 

The 70-100 m/sec zonally symmetric winds in the latitude belt between 30° 
north and 30° south, which are the result of the latitudinal temperature gradient, 
could contribute to the lifting of dust into the atmosphere. However, these winds 
are strongly dependent on the magnitude of the damping parameter employed in 
the model. Although the surface streak patterns observed in the Mariner 9 tele- 
vision pictures (Sagan, et al., 1973) give support to the existance of the predicted 
wind pattern, the actual speed is indeterminate. 

The mechanism producing the diurnal temperature fluctuations is not yet 
fully understood. Tidal resonance amplification does not appear to be a signifi- 
cant factor; except for the winter polar region, dynamic heating is a small part 
of the total atmospheric heating required. While the heating appears to be pri- 
marily a radiative effect, the large infrared atmospheric emissivity required to 
produce adequate cooling at night does not seem to be entirely consistent with 
the observed infrared spectra (Hanel, et al., 1972). 


26 



The linearized dynamic model employed in this study assumes a smooth 
spherical planet and treats frictional forces only through the use of linear 
damping. However, such a model should be adequate for estimating the gross 
properties of the wind fields from the observed temperature fields, and the 
results should provide guidance in the construction of more sophisticated 
models. Future models should include a more realistic treatment of the effects 
of frictional forces on the polar symmetric wind field. Because of the substan- 
tial large scale topographic relief on Mars, orographic effects may play a sig- 
nificant role in atmospheric dynamics (Blumsack, 1971; Gierasch and Sagan, 
1971; Sagan, et al. 1971; Blumsack, et al., 1973), and these effects should also 
be included in future dynamical studies. 
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Table 1. Coefficients of the Legendre functions in the normalized 
geopotential function. Normalization is with respect to RT 0 . The 
a -averaged coefficients are normalized to RT 0 . 


a 

B* 

B* 

^02 

B* 

^01 

B * 

^ n 

1.0 

0.0 

0.0 

0.0 

0.0 

0.9 

-0.0206 

-0.0069 

0.0054 

-0.0024 

0.8 

-0.0434 

-0.0146 

0.0120 

-0.0052 

0.7 

-0.0689 

-0.0233 

0.0199 

-0.0083 

0.6 

-0.0979 

-0.0333 

0.0297 

-0.0120 

0.5 

-0.1314 

-0.0449 

0.0420 

-0.0163 

0.4 

-0.1714 

-0.0589 

0.0577 

-0.0215 

0.3 

-0.2212 

-0.0766 

0.0788 

-0.0282 

0.2 

-0.2879 

-0.1005 

0.1093 

-0.0374 

0.1 

-0.3925 

-0.1384 

0.1607 

-0.0522 


B* 

B * 
h 02 

B* 

CN 

* 


-0.1709 

-0.0595 

0.0638 

-0.0221 
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Table 2. Real and imaginary parts of the coefficients in the Legendre 
function expansion of the normalized pressure ratio. The relationship 
between the diurnal pressure modes (1,1) to (1,6) and the averaged 
geopotential modes (1,1) and (1, 2) are shown with and without damp- 
ing. 




With 

Without 



Damping 

Damping 

Re 

i 

ii 

V 

0.6141 £*j 

- 0.6254 

Re 

* 

a !2=- 

0.6387 /3* 2 

0 

Re 

* 

a !3 = 

0.1054 ^ 

0.1095 fi* n 

Re 

* 

a !4 = 

0.0681 ^* 2 

0 

Re 

A 

Ln 

1 

0.0061 

- 0.0067 /3* x 

Re 

1 

II 

lO 

V 

0.0030 /3* 2 

0 

Im 

a* = 

li 

0.1482 /3*j 

0 

Im 

* — 
a i2~ 

0.4520 ^* 2 

0 

Im 


0.0025 ~P\ X 

0 

Im 


0.0272 /3* 2 

0 

Im 

1 

11 

in 

V 

0.0016 

0 

Im 

A 

II 

0.0004 ^* 2 

0 
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Figure 4. Zonally symmetric mode of the wind field. The product of the ordinate and the mode 
amplitude shown in Figure 7 gives the wind speed. Z Q1 represents the zonal wind 
and M 0] the meridional wind which is the result of damping in the momentum equa- 
tion. 
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Figure 5. Symmetric mode (1, 5) of the zonal wind. The solid line is the temperature phased component. The 
dashed line is the out of phase'component due to damping. The product ofthe ordinate and the mode 
amplitude shown in Figure 8 gives the mode wind speed. 
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Figure 6. Similar to Figure 5 but for the meridional wind. 








METERS/SEC. 

Figure 8. Amplitudes of the diurnal wind field modes. The sign of the real parts of and A^ change between 
cr = .5 and .6 indicating the reversal of the modal wind field. The imaginary parts of all the mode am- 
plitudes, shown dashed, are identically zero if no damping is assumed. 
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Figure 9. Near surface winds during the dust storm. The wind field is generated by the model temperature field 
based on the measured temperatures and the associated surface pressure. Wind amplitudes are scaled 
to the labeled vector. 
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Figure 11. High altitude winds under same conditions as in Figure 9. The wind directions 
are generally reversed from those near the surface. 
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Figure 13. Total heating in an atmospheric column. The maximum heating is at the subsolar point and the total 
heating appears to be direct solar heating. The deviation in the north polar region is due to dynam- 
ics. The heating is expressed in ergs/cm 2 /sec x 10* 5 and degrees K per day in parentheses. 
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Figure 14. Heating due to the vertical motion of the atmosphere. This is in phase with the solar heating and is 
roughly 20% of the total. The heating is expressed in ergs/cm 2 /sec x 10* 5 and degrees K per day 
in parentheses. 
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Figure 15. Heating due to the surface pressure. This is in phase with the solar heating and is roughly 10 o of 
the total heating. The heating is expressed in ergs/cm^/sec x 10*^ and degrees K per day in parentheses. 




